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I.  SUMMARY 


A  fundamental  study  of  the  mechanical  properties  of  microelectronic  thin 
film  materials  has  been  conducted  at  Stanford  University.  The  work  was 
supported  by  AFOSR  Grant  No.  66-0051 .  In  this  final  technical  report  we 
summarize  all  of  the  findings  of  this  research  program  by  citing  the  papers  and 
oral  presentations  that  have  resulted  from  the  work.  The  interested  reader  is 
referred  to  these  published  papers  for  detailed  descriptions  of  this  work.  In  the 
following  sections,  we  describe  some  of  the  work,  done  in  the  last  year,  that  has 
not  yet  appeared  in  print.  Much  of  the  recent  progress  described  in  this  report 
will  serve  as  a  basis  for  the  work  of  the  follow-on  grant. 

It  is  appropriate  to  note  that  before  this  research  program  began,  very  little 
work  had  been  done  in  universities  on  the  mechanical  properties  of 
microelectronic  thin  film  materials.  As  a  result,  much  of  our  early  work  involved 
the  development  of  experimental  techniques,  such  as  sub-micron  indentation 
and  wafer  curvature  measurement,  to  study  stresses  and  mechanical  properties 
of  thin  films  on  substrates.  Interest  in  these  developments  is  indicated  by  the 
accepted  publications  and  invited  oral  presentations  listed  at  the  end  of  this 
report.  The  recent  MRS  Symposium  on  "Thin  Films:  Stresses  and  Mechanical 
Properties"  was  an  outgrowth  of  the  research  work  started  under  this  grant. 

Throughout  this  program  we  have  been  interested  in  the  mechanical 
properties  of  all  kinds  of  microelectronic  thin  film  materials.  We  have  studied 
stresses  and  mechanical  properties  in  interconnect  metals,  passivation  glasses 
and  semiconductors.  Most  of  the  experimental  work  to  date  has  focused  on 
interconnect  metals  and  passivation  glasses.  More  work  will  be  done  on 
semiconductors  in  the  follow-on  grant.  The  wafer  curvature  technique,  that  has 
proved  so  useful  in  the  study  of  metals  and  passivation  glasses,  will  be  used  in 
the  study  of  semiconductor  thin  films.  In  the  present  report  we  present  an  up¬ 
date  on  our  work  on  the  wafer  curvature  technique.  This  includes  a  description 
of  the  factors  that  control  the  sensitivity  of  this  technique  as  well  as  a  report  on 
the  efforts  we  have  made  to  develop  better  facilities  for  making  these 
measurements. 
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A.  Some  Issues  Concerning  Curvature  of  a  Plate  Relevant  To  Thin- 
Film  Stress  Analysis 

(F.  von  Preissig) 

1.  Introduction 

Thin  solid  films  of  many  kinds,  deposited  or  grown  onto  substrates,  are  used 
in  today's  microelectronic  devices  and  for  other  purposes.  Mechanical  sti  esses 
in  these  films  arise  from  various  causes,  including  thermal  expansion 
differences  between  film  and  substrate  materials,  the  physical  nature  of  the 
as-deposited  material,  and  chemical  or  structural  changes  occurring  in  the  films 
upon  processing  after  deposition.  These  stresses  can  lead  to  cracking  or 
delamination  of  the  films  and  can  induce  dislocation  formation  in  the  substrate. 
The  study  of  film  stress  is  also  of  interest  as  a  probe  of  the  physico-chemical 
nature  of  the  film  material. 

In  order  to  experimentally  determine  mechanical  stresses  in  thin  films  on 
substrates,  it  is  common  practice  to  use  the  formula 

~  1  6(1  -  v  )oh 


where 

C  =  curvature  of  substrate;  R  =  radius  of  curvature 
E  =  Young's  modulus  of  substrate 
v  =  Poisson’s  ratio  of  substrate 
h  =  film  thickness 
t  =  substrate  thickness 
cr  =  film  stress. 

Equation  (1)  has  been  derived1-2  for  the  following  case:  The  substrate  is  a  plate 
that  is  thin,  elastically  isotropic,  and  (when  bare)  flat.  The  single  film  has 
uniform  thickness,  much  less  than  f,  and  a  uniform,  isotropic  plane  stress.  The 
temperature  is  uniform.  The  maximum  deflection  due  to  bending  is  less  than 
about  t/2.  The  composite  plate  (substrate  plus  film)  is  mechanically  free. 
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The  films  present  in  integrated  circuit  structures  have  irregular  geometries 
that  complicate  the  understanding  of  mechanical  behavior.  In  particular,  it  is 
necessary  to  have  an  understanding  of  the  stresses  that  exist  in  these 
structures.  Here  we  report  the  results  of  a  comprehensive  study  of  the  stresses 
in  interconnect  lines  bonded  to  substrates.  The  results  are  presented  in  a  form 
that  should  permit  them  to  be  used  for  a  wide  variety  of  thin  film  stress  problems 


Hence,  by  measuring  the  radius  of  curvature  of  a  film/substrate  composite  (e.g., 
by  an  optical  lever  technique),  film  stress  can,  in  principle,  be  calculated. 

In  practice,  however,  conditions  arise  that  completely  invalidate  Eq.  (1), 
require  it  to  be  modified  somewhat,  or  simply  call  its  accuracy  into  question, 
pending  further  analysis.  For  example,  if  the  substrate  is  initially  not  quite  flat,  or 
if  it  is  subjected  to  small,  reproducible  deformations  other  than  that  caused  by 
the  film  stress,  Eq.  (1)  may  apply  only  after  the  non-fiim-stress  effects  are 
subtracted  out,  e.g.,  by  making  a  reference  measurement  of  the  bare-substrate 
surface  profile.3  Also,  thick-film  and  multiple-film  cases  are  covered  by 
extensions  to  Eq.  (1)  that  have  already  been  derived.4-5  In  this  report,  I  examine 
(mostly  theoretically)  some  of  the  effects  of  gravity,  substrate  shape,  film 
nonuniformity,  substrate  crystallinity,  and  temperature  gradients  on  substrate 
curvature. 


2.  Geometric  Conventions,  Definitions,  and  Basic  Formulas 

Geometric  conventions  for  the  composite  plate  as  defined  in  Fig.  1  are 
adopted.  Radial  symmetry  about  the  z  axis  is  assumed  unless  otherwise 
indicated. 


FIG.  1 .  Geometry  (exaggerated,  diametral  cross  section)  of  a 
bent,  circular  plate  with  film,  a  -  radius,  r  -  radial  position,  w  - 
deflection,  C  -  level,  and  <p  -  slope.  Not  shown  are  the  cylindrical 
coordinate  dand  Cartesian  x  and  /  axes. 

The  radial  slope  of  the  deflection  surface  is  given  by 


r>l' 


dw 

dr 


For  small  deflections,  the  local  curvature  is 


d  *w 
dr2 


(2) 


(3) 


(shown  positive  in  Fig.  1). 

If,  in  addition,  C  is  constant  over  r,  then 


and 


<P(r)  -  f  C(r)6r  =  Cr 

o 

1 

=  2a  C- 


(4) 

(5) 


The  substrate,  unless  otherwise  noted,  is  considered  to  be  a  thin  plate 
subjected  to  small  deflections,  as  defined  below  (from  Ref.  6): 

In  the  theory  of  thin  plates,  it  is  customary  to  make  the  following 
assumptions:  (1)  The  plate  is  initially  flat.  (2)  The  material  is 
elastic,  homogeneous,  and  isotropic.  (3)  Thickness  is  small 
compared  to  area  dimensions.  (4)  Slope  of  the  deflection 
surface  is  small  compared  to  unity.  (5)  Deformation  is  such  that 
straight  lines  initially  normal  to  the  middle  surface  remain  straight 
and  normal  to  that  surface.  (Vertical  shear  strains  are  neglected.) 

(6)  Strains  in  the  middle  surface,  arising  from  the  deflection,  are 
neglected  compared  to  strains  due  to  bending.  (7)  Deflection  of 
the  plate  occurs  by  virtue  of  displacement  of  points  in  the  mid¬ 
surface  normal  to  its  initial  plane.  (8)  Direct  stress  normal  to  the 
middle  surface  is  neglected.  (9)  Near  edges  and  boundaries  of 
loaded  areas,  stress  resultants  rather  than  detailed  stress 
patterns  are  considered. 

The  flexural  rigidity  of  the  plate  is  defined  as 

H - 

12(1  -  v2) 
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(6 


The  biaxial  modulus  is 


e=-*-. 

1-v  (7) 

The  state  of  stress  at  a  point  is  called  equi-biaxial  if  there  exists  a  Cartesian 
coordinate  system  for  which  two  of  the  normal  stress  components  are  equal,  the 
third  is  zero,  and  all  shear  stress  components  are  zero;  there  is  simply  an 
isotropic  stress  acting  in  one  plane.  In  this  report,  biaxial  stress  will  refer  to  this 
isotropic  stress  when  it  is  acting  within,  and  uniform  over,  a  given  surface.  The 
film  is  assumed  to  have  biaxial  stress  acting  parallel  to  the  substrate  surface 
and  having  an  average  value  o through  the  film  thickness. 

The  effect  of  the  film  on  the  composite-plate  thickness  and  rigidity  is 
negligible  since  h  «  f  will  always  be  assumed.  The  film  stress  is  positive  when 
tensile,  negative  when  compressive.  Bending  moments  per  unit  length,  M,  are 
positive  if  they  tend  to  cause  positive  curvature  in  their  plane  of  action. 


3.  Effect  of  Gravity 

In  real  situations,  the  sample  is  not  free,  but  is  typically  supported  from 
below  and  subjected  to  the  downward  pull  of  gravity.  The  manner  in  which  this 
uniform  force  affects  the  curvature  depends  on  the  position  of  the  support 
points.  If  the  support  points  are  reproducible,  as  they  are  in  a  tripod 
arrangement  of  pins,  then  the  effect  of  gravity  on  a  given  sample  is  constant  and 
can  be  subtracted  out  by  making  a  bare-substrate  measurement.  If,  however, 
the  sample  sits  on  a  flat  surface  (sometimes  a  tripod  support  is  impractical),  the 
support  points  will  depend  on  the  sample's  warp.  In  this  case,  it  is  desirable  to 
have  an  estimate  of  the  maximum  change  in  curvature  associated  with  a 
change  in  support  points. 

Take  as  two  extremes  a  horizontal,  round  plate  (1)  simply  supported  from 
below  all  around  its  perimeter  and  (2)  supported  only  at  its  center.  With  q 
defined  as  the  uniform  downward  force  per  area,  the  deflection  of  the  plate  for 
the  edge-supported  case  is  given  by6 


(8) 


w. 


g . edge 


5  +  V 

1  +  V 


The  deflection  for  the  center-supported  case  can  be  calculated  by 
superimposing  the  solutions  to  two  different  loading  situations,  as  shown  in 
Figure  2. 


minus 


_1__ 


^yyvi^y/y/  yy 


M 
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/  /  /  / /v  /  /  /  /~y 


FIG.  2.  Demonstration  of  the  way  in  which  the  center-supported 
elastic  plate-bending  case  is  derived  from  the  superposition  of 
two  other  cases. 


The  result  is 


Ly 

0 .  center 


■W 


3+  v 
1+  v 


(9) 


The  difference  between  the  edge  and  center  cases  is 


(10) 


The  corresponding  local  slope  and  curvature,  after  substituting  pgt,  where  p 
is  the  substrate  mass  density  and  g  is  gravitational  acceleration,  are 
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(11) 


(12) 


(The  mathematical  singularity  at  r-  0  makes  Eq.  (12)  physically  unrealistic  very 
near  the  center.) 

The  average  of  A Cg  from  r-  0  to  r=  ron  the  wafer  is 


AC  (r)  = 


j  (a C,(r)  dr 

=  ri HW  -  A0„fr  -»o;  ] 

{,“(,+  ‘'Mi)]- 


3pga: 


er 


(13a) 


For  the  whole  radius  (or  diameter,  due  to  symmetry),  this  average  reduces  to 


AC.  (a)  = 


3  pga‘ 
V  Bt2 


(13b) 


Typically,  however,  what  ci  'ature  measurement  techniques  actually  record  is 
an  average  curvature  that  is  calculated  as  the  slope  of  a  linear  regression  best- 
fit  line  for  a  series  of  measured  substrate-surface  slope  vs.  diametral-position 
data  points.  If  the  curvature  is  not  constant,  then  the  measured  curvature,  Cm, 
will  not  in  general  equal  the  true  average  (as  used  in  Eq.  (13a)]  of  C  from  0  to  r, 
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C  (r).  The  theoretical  value  of  Cm  can  be  derived  from  C(r).  Taking  C(r)  to  be 
the  gravity  effect  ( ACg(r ))  given  by  Eq.  (12)  and  using  ten  slope  measurements 
equally  spaced  on  a  wafer  diameter  to  calculate  the  measured  effect,  ACg.m, 
one  can  derive  that 


ACo  m/  AC0  (a)  =  1.  237  +  0.  237v  , 


(14) 


which  equals  1.27  for  a  quartz  substrate  (v=  .16)  and  has  a  maximum  possible 
value  of  1 .36  (when  v  =  .5).  The  ratio  Cg>J  AC  g (r)  is  closer  to  1  for 
measurements  spanning  less  than  the  full  diameter  (i.e.  r<  a).  We  will 
subsequently  to  refer  to  the  gravity  effect  in  terms  of  ACg(a),  keeping  in  mind 
that  the  magnitude  of  the  effect  as  actually  measured  would  be  comparable. 


In  determining  the  film  stress,  we  may  view  the  maximum  gravity  effect  as  a 
component  of  the  experimental  uncertainty.  With  Co  defined  as  the  curvature 
due  to  film  stress  alone,  as  given  by  Eq.  (1 ),  then 


ACB  (a)  pga2 
C0  -  2ch  • 


(15) 


Note  that  this  relative  variation  is  independent  of  substrate  thickness  and  elastic 
moduli,  but  is  a  strong  function  of  diameter.  To  reduce  the  relative  gravity  effect, 
one  can  simply  use  smaller  substrates. 


Take  the  example  of  a  silicon  (pg  =  2.33  g/cm3  x  980.7  dyne/gram-force) 
substrate  with  a  typical  film  having  o=  200  MPa  *  2x1 09  dyne/cm2  and  h  =  .5 
pm  =  5x1 0-5  cm.  If  the  substrate  is  a  1 00-mm  wafer  (a  =  5  cm),  then  Eq.  (1 5) 
gives  a  relative  uncertainty  of  29%.  But  a  =  1  cm  yields  a  very  tolerable 
uncertainty  of  1.1%.  For  a  substrate  of  quartz  (p  =  2.2  g/cm3,  8  =  8.6x1 011 
dyne/cm2,  v=  .1 6),7  the  relative  gravity  effect  is  slightly  less  than  it  is  for  silicon, 
although  the  absolute  effect  is  larger. 


Experimental  confirmation  of  the  gravity  effect  theory  is  shown  in  Fig.  3.  In 
this  case,  what  is  plotted  is  -d(Awg)/dx  vs.  x,  where  x  is  the  diametral  position, 
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with  the  center  of  the  substrate  being  at  x  =  0.  The  appropriate  theoretical 
equation,  a  modification  of  Eq.  (11),  is 


Experimentally,  slope  vs.  position  data  for  a  round  quartz  substrate  were 
obtained  using  a  laser  reflection  technique.  The  slopes  that  were  measured 
while  the  substrate  was  supported  at  its  center  by  a  2-mm  square  were 
subtracted  (for  each  position)  from  the  slopes  measured  when  the  substrate 
was  supported  around  its  edge  by  a  foam  ring.  The  resulting  experimental 
values  fit  the  theoretical  curve  quite  well. 


Position  (mm) 


FIG.  3.  Experimental  confirmation  of  the  gravity  effect  on  a  quartz 
substrate. 


It  should  be  noted,  however,  that  plate  deflection  is  very  sensitive  to  the 
boundary  conditions,  and  that  in  this  case,  a  uniform,  simple  edge  support  was 


hard  to  achieve.  This  difficulty  resulted  in  some  measurements  varying  as  much 
as  30%  from  those  shown.  In  fact,  in  an  experiment  in  which  the  edge  of  the 
substrate  was  intentionally  supported  at  only  two  points,  on  the  line  of 
measurement,  the  edge-minus-center  measured  curvature  was  2.i7x10-2  nr1, 
or  2.2  times  what  it  was  with  the  uniform  edge  support.  Hence  the  "two 
extremes"  of  support  for  which  the  equations  of  this  section  have  been  derived 
do  not  represent  the  farthest  extremes  of  gravitational  bending  of  an  elastically- 
isotropic,  round  substrate  (to  say  nothing  of  an  anisotropic  or  rectangular  one). 
Nevertheless,  our  analysis  explains  the  general  magnitude  of  the  gravity  effect 
and  shows  how  it  varies  with  substrate  material  and  size.  Such  knowledge  can 
be  useful  in  the  design  of  experiments  and  the  interpretation  of  their  results. 
Unless  a  tripod  support  or  small  substrates  are  used,  the  gravity  effect  is  likely  to 
be  a  significant  source  of  error  in  determinations  of  film  stress  based  on 
substrate  curvature  measurement. 


4.  Effect  of  Shape 

Many  researchers  doing  experiments  that  employ  Eq.  (1)  use  round 
substrates.  But  the  substrate  does  not  need  to  be  round  for  Eq.  (1 )  to  hold,  as 
the  simple  derivation  presented  here  will  demonstrate. 

A  circular  plate  with  a  uniform  bending  moment  per  unit  length  M  applied  at 
its  perimeter  deforms  into  a  spherical  cap  having  a  curvature 

c_  M  _  12  M 

W(1+  v)  Bt 3  (17) 

and  an  internal  bending  moment  isotropically  and  uniformly  equal  to  M.6 
Hence,  a  plate  of  any  shape  will  deform  in  accordance  with  Eq.  (17)  if  Mis 
applied  to  its  edges,  since  the  same  boundary  conditions  are  applied  to  any 
internal  section  of  the  circular  plate. 

In  a  "strength  of  materials"  analysis  (employing  Saint-Venant’s  principle)  of 
the  effect  that  a  thin,  stressed  film  has  on  a  substrate,  a  force  (F)  and  bending 
moment  ( M)  per  unit  length  acting  on  the  perimeter  of  the  substrate  are  found 
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that  make  the  resultant  force  and  moment  at  the  edge  of  the  film  plus  substrate 
equal  zero,  as  illustrated  in  Fig.  4.  The  bending  moment  per  unit  length  is 
•  approximately 


M  =  jtoh  =  M0. 


(18) 


The  force  per  length  F0  =  -crh  acting  at  mid-plane  does  not  affect  the  substrate 
curvature.  Substitution  of  Eq.  (18)  into  Eq.  (17)  yields  Eq.  (1).  Thus,  under 
conditions  in  which  Eq.  (1 )  holds  for  a  round  substrate,  it  also  holds  for  a  sub¬ 
strate  plate  of  any  other  shape. 


h 
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#  FIG.  4.  Model  of  how  a  thin  film  exerts  the  approximate 

equivalent  of  a  bending  moment  per  unit  length  Mo  on  the  edge 
of  a  substrate. 

However,  the  substrate  must  not  be  narrower  than  about  4 1,  or  it  will  cease 

*  to  bend  like  a  plate.8  If  the  width  is  about  equal  to  t  or  smaller,  beam  theory 
applies  and  Ef3/12  replaces  H  (e.g.,  in  Eq.  (17)). 

Although  the  film  effect  on  plate  curvature  is  independent  of  the  shape  of 
®  the  plate,  the  gravity  effect  is  not;  therefore,  in  experiments  in  which  the  gravity 

effect  is  intended  to  be  subtracted  out,  the  substate  shape  shoufd  be  the  same 
for  primary  and  reference  measurements. 
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5.  Effect  of  Nonuniform  Film  Thickness 


Thin  films  are  sometimes  not  of  uniform  thickness  over  the  whole  area  of  the 
substrate.  One  way  to  determine  the  stress  in  such  films  is  to  cut  out  and 
measure  the  curvature  of  a  section  of  the  sample  on  which  the  film  is 
approximately  uniform.  An  alternative  approach  that  involves  using  appropriate 
averages  of  curvature  and  film  thickness  is  applicable  in  some  cases.  This 
method  is  explained  below. 

If,  as  illustrated  in  Fig.  5,  a  radial  bending  moment  per  unit  length  Mis 
applied  to  an  elastically-isotropic  circular  plate  at  r=  r„  then  for  rt.  r„  the  slope 

is8 


Mr 

0(r)-  2  h  ' 


’  1-  v(r.  \21  T  /  r.  \2T) 

}  +  i  +  v  v  a  )  J  L1  v  r  J  Jj 


For  r<  r«, 


Mr  J  i-  v 

«r>  “  Tfo+T)' 1  - 


(19a) 


(19b) 


FIG.  5.  Radial  bending  moment  per  length  M  applied  uniformly  to  a 
circular  plate  at  r=  r* . 

In  the  previous  section,  we  saw  how  a  uniform  thin  film  bends  a  substrate  as 
if  a  force  per  unit  length  oh  were  applied  to  the  surface  of  the  substrate  at  its 
edge,  causing  an  edge-applied  moment  M  =  toh/2.  Similarly,  a  circular  film  of 
thickness  h  but  extending  only  to  r=  r.  on  a  larger  substrate  (a  >  r*)  would  apply 
M  =  toh/2,  but  at  r,  (the  film's  edge)  instead  of  at  the  substrate's  edge.  Now  the 
applied  moment  is  not  the  same  as  the  internal  moment,  even  at  r* ,  since 
reactive  moments  are  introduced.  But  the  applied  moment  is  the  correct 
quantity  to  use  as  Min  Eqs.  (19a)  and  (19b). 


A  film  having  a  radially-symmetric  (with  respect  to  a  round  substrate) 
thickness  distribution  h(r)  can  be  thought  of  as  the  superposition  of  an  infinite 
number  of  disks  of  infinitesimal  (positive  or  negative)  thickness  6h(r).  In 
analogy  with  the  macroscopic  case,  each  of  these  disks  exerts  on  the  substrate 
a  bending  moment  equal  to 


dM(r.)  =  -fdh(rJ, 


(20) 


r *  being  the  location  of  the  disk  edge  and  moment  application.  We  can  imagine 
building  the  film  up  sequentially,  starting  with  a  disk  of  radius  a  and  thickness 
h(a),  then  adding  progressively  smaller  disks  of  thickness 


dh(r.)  =  ^r6r,  (6r.<  0). 


(21) 


The  total  combined  effect  of  these  films  on  the  slope  at  the  edge  of  the 
substrate  is 


<p(a)  =  <t>(r  =  a,  r.~  a,  M  =  tch(a)/2) 


J =0 


+  J  <p(r=  a,  r,  =  r„  M  =d M(r. ))  , 

-  o 


(22) 


where  the  $  components  (here  more  properly  Atf>  and  d0)  are  calculated  using 
Eq.  (19a).  Making  the  indicated  substitutions  yields 


,<7Trh-j;f(4)2^.] 


2H(1 

=  (using  integration  by  parts) 
tea 


2W(1 


(23) 
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But 


l  r.h(r.)drt  =  ~  j  rh(r) dr 
a  o  3  o 

1 


n  a  *  o 


J2nrh(r)dr  =  h (a)  , 


(24) 


where  h  (a)  is  the  average  film  thickness  over  the  whole  area  of  the  substrate. 
The  average  curvature  over  a  full  radial  or  diametral  line  is 


0(a)  =  -1  fc(r)6r  =  ±[t(a) -  HO)] - 


HH 

a 


(25) 


Combining  Eqs.  (23),  (24),  and  (25)  yields 


C(a)  = 


toh(a) 
2H  (1  +  v) 


6(1-  v)oh(a) 
Ef 


(26) 


Equation  (26)  is  the  same  as  the  uniform-film  formula  (Eq.  (1)),  but  with  a  line- 
average  curvature  replacing  C  and  an  area-average  film  thickness  replacing  h. 
Rearranged,  Eq.  (26)  can  be  applied  to  estimate  a  of  films  that  are  roughly 
radially  symmetric  in  thickness  on  a  round,  elastically-isotropic  substrate.  (If  a 
itself  is  a  function  of  r,  ah  (a)  is  used  in  the  equation.) 


Although  the  complete  details  are  not  shown  here,  C  is  not  constant  over  r 
when  the  film  is  not  uniform.  Hence,  as  with  gravity  bending,  0(a),  as  defined 
by  Eq.  (25),  is  not  exactly  equal  to  the  slope  of  a  linear  regression  fit  of  <p  vs.  r 
data,  except  when  <p(a)  and  0(0)  are  the  only  points  used. 

It  also  should  be  noted  that  measuring  curvature  over  less  than  the  full 
diameter  and  ignoring  the  effect  that  the  rest  of  the  substrate  has  on  that 
curvature  can  lead  to  erroneous  calculations  of  stress  when  there  is  film 
nonuniformity.  Suppose,  for  example,  that  a  substrate  of  radius  a  has  a  film  of 
thickness  h  from  r-  0  to  r  =  r.  and  is  bare  ( h  =  0)  everywhere  else.  Then,  with  M 
=  toh/2,  Eq.  (19b)  gives  the  slope  for  r<,  r,  as 
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(27) 


so  the  curvature  for  r<,  r*  is  (using  Eq.  (3)) 

cm  =  60^{i_1^1_(^}=c 

Note  that  although  the  curvature  is  constant,  in  this  case,  over  the  area  where 
the  film  is  (it  varies  for  r>  r»),  its  magnitude  is  smaller  than  that  given  by  Eq.(1) 

by  the  fraction 


Q  = 


1-  v 


(29) 


For  v=  .16  (quartz)  and  rja  =  .5,  this  curvature  reduction  is  31.5%. 


•  FIG.  6  shows  experimental  verification  of  these  equations.  The  bending  of  a 
100-mm  diameter  barium  borosilicate  glass  substrate  (v=  .28)  having  a  50-mm 
diameter  circle  of  silicon  nitride  film  at  its  center  was  compared  to  the  bending  of 
the  substrate  covered  uniformly  by  the  same  film.  The  gravity  effect  and  intrinsic 

•  substrate  curviness  were  subtracted  out  of  the  slope  vs.  x  data  by  means  of  a 
bare-substrate  reference  measurement.  Equation  (26)  predicts  that  the  average 
curvature  of  the  nonuniform  sample  should  be  one-fourth  that  of  the  uniform 
sample,  since  its  h  (a)  is  one-fourth  as  large.  The  experimental  ratio  (Cb/Ca.  for 

•  a  scan  slightly  smaller  than  the  full  diameter)  is  .26.  Equations  (28)  and  (29) 
predict  that  the  curvature  over  the  central  50  mm  of  the  nonuniform  sample 
should  be  constant  and  27%  smaller  than  the  uniform  sample's  curvature.  The 
experimental  value  [(Ca-Cc)/Ca]  is  29%,  with  the  curvature  confirmed  to  be 

•  constant.  Within  experimental  error,  the  experimental  and  theoretical  results 
agree. 
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•  FIG.  6.  Experimental  measurements  of  nonuniform  film  effects. 

Cl  »  curvature  of  the  uniform  sample,  C2  -  average  curvature 
over  96  mm  of  the  nonuniform  sample,  and  C3  -  curvature  of 
central  50  mm  of  the  nonuniform  sample. 

•  Although  Eq.  (26)  was  derived  assuming  radial  symmetry  of  film  and 
substrate,  one  would  also  expect  it  to  apply  to  the  case  of  a  film  that  is  patterned 
in  any  manner  on  a  small  scale  but  uniform  on  a  large  scale  over  a  substrate  of 
any  shape.  However,  when  film  features  are  narrower  than  about  15/7,  there  is 

•  a  significant  amount  of  stress  relaxation  in  the  film,9  so  the  average  value  of  a 
will  not  be  the  same  as  the  pre-patterned  value. 


6.  Effect  of  Crystallinity 

At  the  outset,  I  stated  that  the  derivation  of  Eq.  (1)  includes  the  assumption 
that  the  substrate  is  elastically  isotropic.  Yet  single-crystal  substrates  are  often 
employed  in  thin-film  stress  experiments.  In  this  section,  I  show  that  the 
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anisotropy  associated  with  certain  crystalline  substrates  is  completely 
consistent  with  the  bending  described  by  Eq.  (1 ).  But  for  gravity  bending  and 
nonuniform-film  bending,  the  isotropic-material  equations  are  less  applicable  to 
the  anisotropic  cases. 


The  state  of  stress  in  an  isotropic  plate  with  a  uniform  applied  M  at  its  edges 
(as  discussed  in  Section  4)  is  given  by 


o*(0- 


(30) 


where  om  is  a  biaxial  stress  parallel  to  the  surface  of  the  plate.  When  the 
bending  is  due  to  the  presence  of  a  stressed  thin  film,  then  A/f  =  Mo,  and  there  is 
an  additional  biaxial  stress  <Tfo  =  FqA  (see  Fig.  4).  The  total  biaxial  stress  in  the 
substrate  is  then 


„  (n  -  a  12  f  ,  F°  il 

~  M0  °Fo  ~  t3  ^  t  tit  T 


(31) 


and  there  are  no  other  stresses.  Since  the  slope  of  the  substrate  is  very  small, 
the  stress  state  at  all  points  in  the  substrate  is  described  approximately  by 


CT*)=  o>(s  *,)-  CT,.o 

O3  (=  °;)  =  V)=  °s  (=  T«)=  °6(=  V)  =  °- 


(32) 


Application  of  Hooke's  law  yields  the  corresponding  strain  state: 


°.,o 

€l  =  ^  =  ~B~  ~ 
^..0 

£*~  E 

^4=^=^=0- 


(33) 


The  stress  state  given  by  Eqs.  (31)  and  (32)  is  correct  and  unique  (in  the 
linear  theory  of  elasticity)  because10  (1)  it  satisfies  the  equations  of  static 
equilibrium,  (2)  it  satisfies  the  traction  boundary  conditions  given11  for  all 


surfaces,  and  (3)  the  strains,  related  to  the  stresses  by  the  appropriate 
constitutive  equation  (isotropic-material  Hooke’s  law),  satisfy  the  equations  of 
continuity  (compatibility).12 


The  generalized  form  of  Hooke's  law  for  elastic  anisotropic  media  is13 


£>  =  sfi 


(i  ,j  =  1,2,. ..6), 


(34) 


in  which  the  s,y  (=  s;/)  are  the  elastic  compliances  of  the  material.  In  using  this 
equation,  we  shall  assume  that  the  x  and  y  axes  are  parallel  to  the  plate  at  an 
otherwise  arbitrary  initial  orientation  with  respect  to  elastic  symmetry  elements 
in  the  plate. 


The  state  of  stress  in  an  anisotropic  substrate  bent  by  a  film  is  identical  to 
that  given  [Eq.  (31 )]  for  an  isotropic  substrate  because  this  state  again  satisfies 
the  conditions  required  for  a  unique  solution.  The  boundary  conditions  and 
equilibrium  equations  are  the  same  as  in  the  isotropic  case,  so  they  are 
satisfied.  The  equations  of  continuity,  which  equate  sums  of  second  derivatives 
of  strain  with  respect  to  position,  are  satisfied  because  the  strain  components 
are  simply  linear  combinations  [by  Eq.  (34)]  of  the  stress  components,  which 
vary  linearly  with  position14  (making  all  second  derivatives  zero). 


The  state  of  strain  in  the  anisotropic  plate  remains  to  be  determined. 
Substituting  the  stresses  from  Eq.  (32)  into  Eq.  (34)  yields,  for  /'  =  1  and  2, 


*1=(S11  +  S12)°,.0 
£Z  =  (®12  +  ®22)°*.CT 


(35) 


Converting  the  compliances  to  engineering  elastic  constants,15 

S,,+  s,2=  5i(1  +  ^)  =  -^(1“ 

s12+  s22=  S22(^1  +  — j  =  -^(1-  v21)  =  —  . 


(36) 
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The  rest  of  the  strain  components  are 


^3  ~  (^13  +  ®23)°»,0 
*4=(S14  +  SmKo 

^=(*15  +  S2s)°..0 
^  =  (  ®16  +  0  ' 


(37) 


For  materials  that  have  third  order  or  higher  rotational  elastic  symmetry 
("polysymmetry")  about  the  zaxis,  the  strain  state  is  greatly  simplified.  For  these 
materials,  S22  =  sn-16  Thus,  from  Eqs.  (35)  and  (36), 


£,=  £?  = 


(38) 


Also,  S23  =  S13,  so 


2v. 


11 


(39) 


Finally,  the  compliances  are  such  that 

e4=^=^  =  °-  (40) 

This  strain  state  is  precisely  analogous  to  that  given  by  Eq.  (33)  for  an  isotropic 
substrate,  with  By 2  replacing  B  and  S13  replacing  -v/E.  ey  and  cts,o  are  constant 
for  rotations  of  the  coordinate  system  about  the  z  axis,  hence  Eq.  (38)  implies 
that  812  is  similarly  invariant.  (S13  is  invariant  for  z-axis  rotations  for  all 
materials.)  Thus,  for  substrates  having  z-axis  polysymmetry,  the  curvature  of 
any  constant-^  surface  is  uniform  and  given  by  an  analogue  of  Eq.  (1), 
specifically 


where  812  is  the  in-plane  biaxial  modulus  of  the  plate.  Values  of  S12  for  four 
cubic  semiconductor  crystals  in  (100)  and  (111)  wafer  orientations  are  given  in 
Ref.  17. 

The  deflection  of  a  simply-supported  anisotropic  circular  plate  subjected  to 
loading  normal  to  its  surface  is  extremely  difficult  to  solve  analytically.1®  Hence, 
I  do  not  know  the  precise  effect  of  crystal  anisotropy  on  gravity  bending. 
Certainly,  Eq.  (12)  cannot  describe  the  deflection  of  (100)-Si  (fourth-order 
symmetric)  since  the  equation  contains  v  separate  from  8,  and  the  in-plane  v 
(vi2)  varies  with  6  for  this  substrate.17  Even  for  (1 1 1  )-Si  (third-order  symmetric), 
which  has  isotropic  in-plane  E  and  v  (also  Sn  and  s  12),  a  radial ly-sym metric 
deflection  is  doubtful  because  shear  stresses  of  the  type  05  exist  and  S15  varies 
with  6.  Nevertheless,  for  the  purpose  of  estimating  the  magnitude  of 
experimental  uncertainty  due  to  the  gravity  effect,  the  equations  of  Section  3 
(perhaps  using  directionally-averaged  elastic  moduli)  should  be  adequate  for 
most  substrate  materials. 

The  nonuniform-film  deflection  equations  presented  in  Section  5  also  fail  to 
take  into  account  the  anisotropy  of  (100)-Si  v12.  But  they  are  probably  nearly 
exact  for  (1 1 1)-Si,  since,  unlike  the  case  of  gravity  loading  mentioned  above, 
only  01  and  <r2  are  nonzero  (at  least  for  an  isotropic  substrate). 


7.  Effect  of  Thermal  Gradient 

In  some  experiments,  a  film  on  a  substrate  is  heated  while  the  curvature  is 
being  measured  in  order  to  determine  the  effect  of  temperat  ire  changes  and 
annealing  on  film  stress.  Thermal  gradients  may  exist  in  the  substrate,  and  their 
effect  on  the  curvature  must  be  considered. 

We  will  examine  the  case  in  which  the  substrate  has  a  temperature  that 
varies  linearly  from  7  at  the  top  surface  to  7+  A  Tat  the  bottom  surface.  This 
vertical  thermal  gradient  contributes  directly  to  curvature,  since  the  thermal 
strain  varies  linearly  with  £.  The  thermally-induced  curvature,  which  is  uniform 


r  -  «  AT 

'"/th  f  > 


where  a  is  the  thermal  expansion  coefficient  of  the  substrate. 

In  the  following  examples,  we  will  use  the  material  parameters  given  in 
Table  1.19.20 


Table,  1 


Thermal  Exp.  Coef.  (10-®  °C*1)  Thermal  Conductivity  (W/m°C) 


a,  25  °C 

a,  700  °C 

A,  300  °C 

A,  700  °C 

Silicon 

2.6 

4.4 

136 

32 

Fused  Quartz 

.50 

.39 

1.7 

2.7 

Suppose  we  wish  to  know  what  AT  would  cause  a  Cth  of  IxlO  3  nr1,  which 
is  a  small  but  significant  curvature  change,  in  a  .5  mm-thick  substrate.  From  Eq. 
(42),  AT  =  Cihfcfa  = 


for  silicon: 
for  quartz: 


.19  C°  at  25  °C  ,  .114  C°  at  700  °C; 
1 .0  C°  at  25  °C  ,  1 .28  C°  at  700  °C. 


Let  us  calculate  how  much  power  must  pass  vertically  through  a  round  100- 
mm  wafer  to  cause  the  above  effects. 

Heat  flow  rate  =  IA~-  ; 

Ax  (43) 

A  »  k(.05  m)2,  Ax  =  t  =  5x1 0-4  m.  The  required  power  is 

for  silicon:  406  watts  at  25  °C,  57  watts  at  700  °C; 


for  quartz: 


27  watts  at  25  °C 


54  watts  at  700  0  C 


Hence,  for  silicon  substrates  at  high  temperatures  and  quartz  substrates  at  all 
temperatures,  vertical  thermal  gradients  within  the  substrate  will  have  a 
significant  effect  on  curvature  unless  the  furnace  used  has  a  low  vertical  heat 
flux  within  the  compartment  holding  the  substrate. 

Horizontal  thermal  gradients  may  also  exist  in  the  substrate.  They  induce 
thermal  stresses  that  can  alter  the  curvature  of  an  already-curved  substrate,  but 
the  analysis  of  this  effect  will  not  be  attempted  here. 


8.  Summary 

We  have  seen  that  the  standard  equation  relating  substrate  curvature  to  film 
stress  is  theoretically  accurate  under  certain  experimental  conditions  and  not 
accurate  under  others.  Gravitational  forces  can  cause  a  large  error,  but  one  can 
circumvent  this  problem  by  using  small  substrates  or  reproducible  supports. 

The  shape  of  the  substrate,  as  long  as  it  is  a  thin  plate,  does  not  affect  the 
curvature-stress  relation.  Films  of  nonuniform  thickness  are  amenable  to 
analysis  by  a  simple  analogue  of  the  classical  curvature-stress  relation  in 
certain  cases.  The  crystallinity  of  (1 1 1 )-  or  (1 00)-oriented  cubic-material 
substrates,  among  others,  necessitates  another  analogue.  Thermal  gradients 
can  significantly  alter  substrate  curvature.  These  pieces  of  information  should 
be  valuable  to  researchers  measuring  stresses  in  thin  films  on  substrates. 
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B.  Redesigned  Furnace  and  Ambient  Gas  Control  for  Thin-Film 
Stress  Measurement  System  (F.  von  Preissig) 


1.  Background 

Currently,  we  have  an  apparatus  that  measures  mechanical  stress  in 
thin-film  materials  on  substrates.  Film  stress  is  calculated  from  optically- 
measured  substrate  curvature.  Ours  is  one  of  the  very  few  systems  that  attempt 
to  heat  the  sample  while  measuring  its  stress.  This  permits  in-situ  monitoring  of 
thermally-induced  changes  in  the  film  material  that  are  reflected  in  the  changing 
stress.  Although  our  system  has  proved  to  be  an  important  research  tool,  its 
heating  chamber  has  certain  inadequacies  and  limitations  that  make  it 
unsuitable  for  the  new  generation  of  experiments  that  should  be  done.  I  have 
designed  and  am  building  a  new  furnace  that  should  eliminate  the  problems 
and  add  new  capabilities.  This  furnace  is  now  near  completion.  In  addition,  I 
have  created  equipment  to  introduce  argon-diluted  steam  into  the  sample 
compartment. 

2.  Problems  With  the  Current  Furnace 

The  sample  substrate  sits  directly  on  a  flat  graphite  susceptor.  Since  the 
substrate  is  not  perfectly  flat,  it  contacts  the  surface  discontinuous^  and  in  an 
unpredictable  way,  causing  two  main  problems.  First,  the  effect  of  gravity  on 
curvature,  which  can  be  large,  is  not  constant,  so  it  cannot  be  subtracted  out 
using  a  reference  measurement.  Second,  the  sample  may  shift  and  wobble 
during  the  experiment,  adding  scatter  to  the  data,  due  to  unstable  seating. 

Since  heat  is  supplied  only  from  below,  large  thermal  gradients  exist  in  the 
sample  compartment.  These  temperature  differences  probably  don’t  affect  the 
sample  temperature  much  when  the  substrate  sits  on  the  graphite  (containing 
the  measurement  thermocouple).  However,  in  order  to  counter  the  gravity  and 
stability  problems,  the  substrate  would  need  to  sit  on  a  tripod  slightly  above  the 
graphite  surface.  But  when  this  is  done,  the  sample  temperature  differs 
significantly  from  that  measured  by  the  thermocouple  (40  C°  at  400  °C).  For 


most  studies  of  kinetic  processes,  even  ten-degree  errors  are  not  tolerable. 
Temperature  gradients  are  also  responsible  for  gas-density  fluctuations,  which 
deflect  the  laser  beam,  and  can  induce  thermal-strain  curvatures  in  the 
substrate. 

The  sample  is  exposed  to  insulation  material  within  the  gas-tight 
compartment.  This  configuration  reduces  ambient  gas  (e.g.,  argon)  purity. 


3.  Description  and  Capabilities  of  the  New  Design 

The  new  furnace,  shown  schematically  in  Fig.  1,  has  a  fused-quartz  inner 
chamber,  including  a  gas-tight  seal.  Hence,  the  ambient  gas  remains  pure,  and 
reactive  gases  such  as  oxygen  and  steam  may  be  introduced. 

Surrounding  the  quartz  on  all  sides  (leaving  an  optical  opening)  is  a 
thermally-conductive  graphite  sheath,  heated  from  above  and  below.  This 
arrangement  will  essentially  eliminate  thermal  gradients  in  the  sample 
compartment. 

The  sample  substrates  rest  on  quartz  tripods  within  the  chamber,  thus 
sustaining  a  constant,  subtractable  gravity  effect  and  attaining  seating  stability. 
Since  there  are  no  thermal  gradients,  the  sample  temperature  is  accurately 
measured  (at  least  in  the  steady  state)  by  a  thermocouple  probe  in  the  sample 
compartment. 

The  furnace  opens  in  such  a  way  that  samples  can  be  inserted  and 
removed  (by  means  of  a  vacuum  chuck  assembly)  quickly  while  the  furnace  is 
hot  (see  Fig.  2).  Samples  will  heat  up  to  a  constant  annealing  temperature  in 
about  a  minute  instead  of  the  10-20  minutes  currently  possible.  Annealing 
temperatures  up  to  990  °C  will  be  possible. 

The  ambient-control  system  shown  and  described  in  Fig.  3  complements 
the  new  furnace's  capability  of  handling  steam  for  use  in  solid-gas  reaction 
experiments. 


FIG.  2.  Furnace  ready  to  accept  sample  carried  by  vacuum  chuck  assembly. 


FIG.  3(a)  Ambient  gas  control  system. 
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FIG.  3(b).  Gas  control  symbols. 


All-Dry  Mode 

For  experiments  requiring  an  inert-gas  ambient,  helium  and  argon  can  be 
used.  Typically,  the  furnace  is  purged  of  air  by  helium  (A)  immediately  after 
sample  insertion.  The  high  thermal  conductivity  of  helium  minimizes  the  time  it 
takes  for  the  sample  to  reach  furnace  temperature.  After  the  temperature  has 
stabilized,  valve  B  is  used  to  switch  from  helium  to  argon  (C)  flow,  and  argon  is 
used  for  the  remainder  of  the  experiment.  The  inert-gas  flow  rate  is  controlled 
by  flowmeter  D.  After  leaving  the  furnace,  inert  gases  percolate  through 
silicone  oil  in  the  bubbler  to  eliminate  back-diffusion  of  air  into  the  furnace.  The 
gas  pressure  in  the  furnace  is  always  approximately  one  atmosphere. 

Steam/Dry  Mode 

Pure  steam  is  produced  when  de-ionized  water  from  a  storage  tank  situated 
above  the  apparatus  flows  through  flowmeter  E  at  a  controlled  rate  and  is 
vaporized  at  the  same  rate  in  the  boiler.  The  steam  can  be  diluted  a  desired 
amount  by  argon  controlled  by  flowmeter  F.  A  separate  argon  cylinder  (G)  is 
used  for  this  purpose  to  prevent  possible  contamination  of  house  argon  by 
water  vapor.  To  reduce  contamination  of  this  cylinder,  a  plug  valve  (H)  is 
closed  when  no  diluent  argon  is  flowing.  Diluent  argon  and  pure  steam  mix  at 
junction  I.  The  portion  of  this  mixed  gas  to  be  used  in  the  long  term  is  controlled 
by  needle  valve  J.  The  rest  either  exits  to  the  exhaust  vent  or  is  shunted  toward 
the  furnace  for  use  as  an  initial  purge,  depending  on  the  setting  of  valve  K.  The 
four-way  valve  (L)  allows  an  instantaneous  switch  to  be  made  between  wet  and 
dry  gases  flowing  through  the  furnace  (with  the  gas  not  currently  selected  going 
directly  to  the  exhaust).  Valve  M  is  used  to  send  dry  gas  through  the  bubbler 
and  wet  gas  directly  to  the  exhaust. 


FIG.  3(c).  Gas  control  operation. 


4.  Planned  Experiments 


Studies  have  begun  on  the  following  material  processes:  (1 )  Densification 
and  flow  in  borophosphosilicate  glass  (BPSG)  films  upon  annealing.  (2)  Water 
diffusion  in  BPSG  films  exposed  to  wet  and  dry  ambients  at  various 
temperatures.  (3)  Crystallization  in  amorphous  silicon  films  upon  annealing.  It 
may  also  be  possible  to  monitor  wet  oxidation  of  silicon  in-situ.  Preliminary 
results  indicate  that  the  associated  stress  changes  are  large  and  can  be  used  to 
learn  about  the  kinetics  of  these  processes. 


C.  Mechanical  Properties  of  Thin  Aluminum  Metallizations 

(J.F.  Turio) 

As  the  dimensions  of  integrated  circuits  continue  to  shrink,  an 
understanding  of  the  mechanical  properties  of  thin  films  becomes  ever  more 
important  to  device  manufacture,  operation  and  reliability.  Several  methods 
have  been  used  to  study  thin  film  stresses.  The  variation  of  film  stress  with  film 
thickness  has  been  monitored  by  measuring  the  deflection  of  thin  cantilever 
beams  which  act  as  the  substrate  during  deposition.  X  -  rays  have  frequently 
been  used  to  measure  the  elastic  strain  of  the  substrate  and  this  can  be  related 
to  the  stress  in  the  film.  Also,  lasers  have  been  used  in  several  ways  to 
measure  sample  curvature,  which  can  also  be  related  to  the  stress  in  the  film. 

At  Stanford,  a  laser  based  curvature  measurement  system  has  been 
designed  for  quick  and  convenient  acquisition  of  thin  film  stress  data.  A  laser 
beam  is  aimed  at  a  mirror  that  is  mounted  on  a  shaft  whose  rotation  angle  can 
be  precisely  controlled.  The  beam  then  travels  to  the  sample.  As  the  shaft 
rotates,  the  laser  beam  sweeps  across  the  sample.  In  addition,  the  beam 
passes  through  a  lens  so  that  the  beam  will  strike  the  sample  at  the  same 
incident  angle  regardless  of  shaft  angle.  After  reflection  from  the  sample,  the 
beam  passes  back  through  the  lens  which  focuses  it  onto  a  position  sensing 
photodiode.  The  sample  curvature  is  obtained  from  a  plot  of  the  shaft  angle 
versus  the  photodiode  signal.  A  single  curvature  measurement  at  room 
temperature  can  be  made  in  about  10  seconds.  The  system  can  also  heat  a 
sample  to  450  °C  at  0  to  10  degrees  per  minute. 

The  system  has  a  theoretical  resolution  of  25x1  O’6  nr1,  which  corresponds 
to  a  surface  with  a  40  kilometer  radius  of  curvature.  In  practice,  the  system 
displays  a  resolution  of  IxlO*3  nr1  or  a  1000  meter  radius  of  curvature.  The 
corresponding  stress  resolution  is  about  10  MPa  but  is  sample  dependent.  The 
overall  system  accuracy  is  roughly  5  %.  A  typical  curve  of  stress  versus 
temperature  obtained  with  this  system  is  shown  in  Figure  1 . 
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Figure  1 .  A  typical  plot  of  stress  versus  temperature  for  a  1  micron 
Al  - 1%  Si  film  on  a  silicon  substrate.  The  open  diamonds 
are  heating  data  and  the  filled  diamonds  are  cooling  data. 


During  heating,  the  stress  changes  by  elastic  deformation  until  a 
temperature  and  compressive  stress  are  reached  that  cause  significant  plastic 
flow.  The  maximum  compressive  stress  attained  is  about  120  MPa  which 
corresponds  to  about  17  ksi.  Bulk  aluminum  samples  typically  yield  between  5 
and  15  ksi.  During  cooling,  no  notable  elastic  range  is  discernible  but,  the  flow 
stress  in  tension  increases  with  decreasing  temperature.  A  change  in  slope  is 
seen  at  about  250°C,  which  could  be  due  to  precipitation  of  silicon.  The  feature 
is  more  apparent  in  samples  that  contain  copper. 

The  results  in  Figure  1  are  complicated  by  the  fact  that  stress  is  a  function 
of  both  temperature  and  time.  But,  isothermal  tests  can  also  be  performed.  A 
typical  stress  relaxation  curve  is  shown  in  Figure  2. 


29 


Figure  2.  Stress  relaxation  data  for  a  1  micron  aluminum  film  on  a 
silicon  substrate. 


The  advantage  of  an  isothermal  test  is  that  the  data  can  be  cast  into  a 
standard  rate  form  and  compared  to  bulk  data.  In  doing  so  for  several  samples, 
the  stress  at  a  given  strain  rate  is  seen  to  be  higher  for  thin  film  aluminum, 
implying  a  higher  resistance  to  plastic  flow  than  in  bulk  samples.  However,  the 
data  is  sparse  and  somewhat  inconsistent.  In  fact,  the  log  of  stress  was  plotted 
against  the  log  of  time  for  many  isothermal  tests  to  determine  the  stress 
exponent  in  the  simple  equation, 

£p“On 

but  the  resulting  exponent,  n,  varied  between  7  and  40. 

The  inconsistencies  in  the  above  data  could  be  caused  by  the  furnace 
which  suffers  severe  wear  and  distortion  during  thermal  cycling.  At  times,  the 
temperature  has  been  found  to  be  inaccurate  by  as  much  as  10  %.  Another  part 
of  the  problem  could  be  that  the  system  can  only  accurately  resolve  a  stress  of 
10  MPa.  The  stress  in  Figure  2  only  varies  80  MPa  in  24  hours  at  one  third  of 


the  melting  point  of  aluminum  while  tests  at  higher  temperatures  may  begin  at 
only  50  MPa  and  vary  by  only  20  MPa.  Lower  temperature  tests  cannot  be 
performed. 

Consequently,  a  system  has  been  designed  and  is  being  built  that  will 
allow  a  wider  range  of  testing  conditions  with  improved  accuracy.  The  design 
specifications  include  heating  and  cooling  capabilities  between  -50°C  and 
800°C  at  -20°C  to  20°C  per  minute.  The  theoretical  resolution  is  over  100 
kilometers  but,  the  actual  resolution  will  be  much  less.  The  expected  stress 
resolution  is  less  than  1  MPa  with  an  accuracy  of  better  than  1  %. 

Typical  factors  influencing  the  mechanical  behavior  of  metals  include 
composition,  grain  size  and  structure  but,  thin  films  have  also  shown 
dependence  on  film  thickness,  presence  of  oxides  and  interface  quality.  The 
stress  in  a  0.2  micron  thick  aluminum  film  has  been  found  to  be  600  MPa,  which 
is  five  times  the  stress  in  a  1  micron  thick  film.  And,  as  noted,  precipitation 
hardening  seems  to  play  some  role  in  film  strength  but  the  effect  has  not  been 
quantified. 

In  an  effort  to  obtain  a  greater  understanding  of  the  mechanisms  governing 
the  strength  properties  of  aluminum  thin  films,  a  more  complete  study  will  be 
conducted  using  the  new  curvature  measurement  device.  General  thermal 
cycles,  as  in  Figure  1,  and  isothermal  tests,  as  in  Figure  2,  will  be  performed 
over  a  range  of  film  parameters  to  determine  the  magnitude  of  their  effect.  In 
addition,  isothermal  results  can  be  related  to  bulk  data  and  established 
deformation  mechanisms. 


D.  Finite  Element  Calculations  of  Thermal  Stresses  in  an 

Unpassivated  Line  Bonded  to  a  Rigid  Substrate 

(A.I.  Sauter) 

1.  Introduction 

Stresses  develop  in  integrated  circuit  structures  during  processing  and 
service  as  a  result  of  a  variety  of  effects.  Thermal  expansion  coefficients  of 
neighboring  materials  may  be  sufficiently  different  that  temperature  changes 
produce  significant  stresses.  Also,  adjoining  materials  may  react  in  such  a  way 
that  the  product  phases  have  a  volume  which  differs  from  that  of  the  reacting 
phases.  This,  too,  results  in  the  generation  of  stresses.  In  addition,  stresses 
develop  naturally  in  epitaxial  films  when  their  lattice  parameters  differ  from  that 
of  the  substrate.  In  short,  there  are  many  possible  sources  of  stresses  in 
integrated  circuits.  These  stresses  can  lead  to  plastic  flow,  cracking,  void 
formation,  and  eventually  failure  of  the  circuit.  It  is  therefore  of  interest  to  know 
the  magnitudes  and  concentrations  of  such  stresses. 

Considerable  attention  has  been  focussed  on  the  stresses  in  continuous 
thin  film  layers.  Stresses  in  thin  films  arising  from  epitaxial,  thermal  and  growth 
mechanisms  have  been  well  studied  and  are  generally  understood.  However, 
thin  film  materials  used  in  integrated  circuit  technology  are  almost  exclusively 
found  not  as  continuous  films  but  as  patterned  films.  Thus  it  is  of  interest  to 
study  the  stresses  as  they  develop  in  irregular,  patterned  structures. 
Interconnect  metals  are  particularly  susceptible  to  thermal  stresses,  because 
their  thermal  expansion  coefficients  are  so  much  greater  than  all  of  the  other 
materials  involved.  The  other  materials  are  usually  covalently  bonded 
insulators  and  semiconductors  and  have  similar  thermal  expansion  coefficients, 
while  the  thermal  expansion  coefficients  of  the  metals  can  be  up  to  an  order  of 
magnitude  larger.  For  example,  one  of  the  most  widely  used  metallization 
materials  is  aluminum.  Its  thermal  expansion  coefficient  is  23.6x1 0*6  /  °C,  while 
that  for  Si  is  about  4x10*6  /  °C. 

It  is  important  to  understand  how  stresses  in  patterned  films  differ  from 
those  in  continuous  films.  As  yet,  the  stresses  in  patterned  films  have  not  been 
measured  or  characterized  in  a  comprehensive  way.  It  is  often  difficult  to 


calculate  the  expected  stresses  under  a  given  set  of  conditions  because  of  the 
often  complicated  geometries  involved.  However,  the  finite  element  method 
can  be  used  to  find  the  thermal  stresses  in  these  cases  [1]. 

2.  The  Finite  Element  Method 

The  finite  element  method  is  used  in  the  present  work  to  find  the 
distribution  of  stress  in  a  long  line  of  film  bonded  to  a  rigid  substrate  and 
subjected  to  a  thermal  strain.  The  geometry  involved  and  a  sample  mesh  are 
shown  in  Figure  1 .  The  thickness  of  the  line  is  h  and  its  width  is  w.  Calculations 
have  been  made  for  width  to  thickness  ratios  of  w/h  =  2,  4,  6,  8,  and  10.  Note 
that  the  y-z  plane  is  a  plane  of  symmetry  for  the  line,  so  the  mesh  only  covers 
the  right  half  of  the  line.  The  left  half  of  the  line  is  simulated  by  the  boundary 
condition  that  there  be  no  x  displacements  on  the  plane  x  =  0.  Also,  two 
dimensional  plane  strain  elements  are  used  in  the  mesh.  The  strain  in  the  z 
direction  is  taken  to  be  zero;  this  naturally  causes  stresses  to  develop  in  the 

line. 


The  MARC  finite  element  program  was  used  to  solve  the  problem  of 
thermal  stresses  in  a  long  thin  line  bonded  to  a  rigid  substrate.  The  meshes 
used  were  generated  using  the  MENTAT  program.  The  material  properties  of  Al 
and  a  temperature  change  of  -475°C  (from  500°C  to  25°C)  were  used  in  the 
calculation,  but  the  results  presented  have  been  normalized  to  make  them 
material  and  temperature  difference  independent.  In  order  to  test  whether  the 
meshes  were  fine  enough  to  give  accurate  results,  the  element  size  was  refined 
by  a  factor  of  two;  the  results  were  not  significantly  different.  It  was  concluded 
that  the  original  meshes  were  sufficiently  fine  to  show  the  important  features  of 
the  stress  distributions  in  the  lines.  All  calculations  were  made  using  a 
VAXstationll/GPX  computer. 


rigid  bond 


Substrate 


ujfO 


y 

L: 
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u^u^O 


Figure  1 . 


(a)  A  schematic  of  the  geometry  of  the  problem.  The  y-z 
plane  is  a  plane  of  symmetry,  (b)  A  sample  mesh  for  an 
aspect  ratio  of  w/h  =  8.  The  boundary  conditions  are 
shown.  Note  that  the  left  half  of  the  line  is  omitted  due  to 
symmetry. 


3.  Plane  Stress  vs.  Plane  Strain 

In  1949,  Aleck  [2]  numerically  solved  the  problem  of  thermal  stresses 
induced  in  a  thin  rectangular  plate  clamped  along  one  edge.  The  current 
problem  is  very  similar  to  Aleck’s  problem,  except  that  Aleck  used  plane  stress 
conditions  for  his  plate  while  the  more  suitable  condition  for  a  long  thin  line  is 
plane  strain.  Plane  stress  calculations  were  also  made  in  order  to  compare 
them  with  the  previous  numerical  results.  Figure  2  shows  the  distribution  of 
transverse  thermal  stress,  oxx,  in  a  line  with  an  aspect  ratio  of  w/h  =  10.  The 

plane  stress  numerical  results  of  Blech  and  Levi  [3]  obtained  using  Aleck's 
method  are  compared  with  the  finite  element  results.  Good  agreement  between 
the  two  methods  is  found.  The  main  difference  between  the  two  is  that  the  finite 
element  results  do  not  show  the  stress  at  the  interface  (y/h  =  0)  returning  to  zero 
at  the  edge  of  the  line,  but  rather  indicate  a  sharp  increase  in  stress.  In  fact  the 
stress  calculated  at  that  point  should  depend  on  the  path  used  to  approach  the 
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point  [4].  It  should  be  zero,  as  Aleck  shows,  if  the  point  is  approached  from 
outside  the  line,  or  along  the  outer  edge  of  the  line.  However,  it  may  have  some 
other  value  if  approached  from  inside  the  line,  as  a  singularity  exists  there.  So 
the  discrepancy  is  probably  due  to  a  difference  in  calculation  paths.  Thus, 
except  for  this  point,  we  conclude  that  the  finite  element  method  can  give  an 
accurate  solution  to  this  problem. 


Figure  2.  A  comparison  of  numerical  and  FEM  results  for  the  plane 
stress  case,  showing  the  transverse  stress  in  a  line  of  w/h 
=  10  aspect  ratio.  The  numerical  results  are  those  of 
Blech  and  Levi  [3]. 


For  both  the  plane  strain  and  plane  stress  cases,  the  stresses  are 
normalized  by  the  value  of  aXx  that  would  exist  in  the  center  of  the  line  at  the 

interface  between  the  line  and  the  substrate  if  the  line  were  very  wide.  These 
values  are  computed  as  follows:  For  the  case  of  plane  stress  (ojz  =  °)<  °yy  =  0 

because  there  is  no  constraint  in  the  y  direction  (Figure  3(a)).  Since  the  line  is 
perfectly  bonded  to  the  substrate, at  the  interface  the  total  strain  in  the  x  direction 
must  be  zero,  hence:  eXXi elastic  =  'exx,thermal-  By  Hooke's  law, 


Lr^  _ %i/ st  \ 

1  ^xx 

£  — 
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where  E  is  Young's  modulus  of  the  film,  and  v  is  Poisson’s  ratio.  The  thermal 
strain  is  given  by  -AaAT,  where  Aa  is  the  difference  in  thermal  expansion 
coefficients  between  the  film  and  the  substrate,  and  AT  is  the  difference 
between  the  initial  and  final  temperatures.  Thus  the  transverse  stress'at  the 
center  of  a  wide  line  is  given  by 


o  xx  =  E  Aa  AT 


(2) 


For  the  case  of  plain  strain,  Oyy  =  0  as  before,  but  Gjz  -  axx*  since  the  middle  of 
a  long  wide  line  should  behave  exactly  like  a  thin  film,  and  therefore  have  a 
biaxial  stress  state  (Figure  3(b)).  As  before, 

lr  -i  (1_v)axx 

8  xx,  elastic  ~  EL°XX  ”  v(ayy  +  azz)j  “  E  (3) 


so  for  the  case  of  plane  strain,  the  expected  stress  is  given  by 


o  =  —  E— rAaAT 
xx  (1  —  v) 


(4) 


Thus,  it  is  expected  that  near  the  center  of  a  wide  line  the  transverse  stress 
for  plane  strain  will  be  larger  than  for  plane  stress  by  a  factor  of  1/(1-v).  This 
prediction  is  borne  out  by  the  finite  element  results,  as  for  both  plane  strain  and 
plane  stress  the  normalized  transverse  stress  approaches  1.0  at  the  center  of 
the  line  as  the  line  width  increases.  In  all  of  the  figures  showing  normalized 
stress,  the  normalization  factor  is  given  by  equation  (2)  for  plane  stress  and 
equation  (4)  for  plane  strain. 


Aleck's  results  have  been  used  by  others  [5,6]  to  describe  stresses  in  thin 
film  structures,  however  these  authors  left  the  calculation  in  plane  stress  form. 
Although  it  is  possible  to  approximate  the  plane  strain  results  by  dividing  the 
plane  stress  results  by  (1-v),  the  exact  form  of  the  stresses  cannot  be  obtained 
by  this  procedure. 


Figure  3.  (a)  Aleck's  geometry,  suitable  for  plane  stress 

calculations,  (b)  The  geometry  of  a  line  of  metallization, 
suitable  for  plane  strain  calculations. 


4.  Results 


The  transverse,  shear,  normal,  and  longitudinal  stresses  (aXx.  <*xy,  ayy  ancl 
azz)  are  shown  for  a  line  with  aspect  ratio  w/h  =  6  in  Figure  4.  They  are  plotted 
for  five  different  values  of  y/h.  With  the  mesh  that  was  used,  it  was  not  possible 
to  calculate  the  stresses  at  y/h  =  0  and  y/h=1 ,  so  the  results  are  given  for  y/h  = 
0.00704  and  y/h  =  0.993  instead.  Thus  the  stresses  at  the  interface  and  at  the 
top  surface  of  the  line  are  not  known  exactly.  The  transverse  stress,  Figure  4(a), 
approaches  its  expected  value  at  the  center  of  the  line  near  the  interface.  The 
stress  drops  to  zero  at  the  edge  of  the  line  (x/h  =  3)  except  near  the  interface, 
where  there  is  a  singularity.The  shear  stress,  Figure  4(c),  is  zero  at  the  top  of  the 
line,  becoming  larger  as  one  approaches  the  interface  and  having  a  maximum 
at  each  y  value  near  the  edge.  At  the  interface,  the  shear  stress  appears  to 
have  a  singularity  at  the  edge.  Again,  this  may  be  a  matter  of  calculation  path, 
or  it  may  be  that  there  is  a  turn  toward  zero  in  this  curve  that  is  so  near  the  edge 
it  is  impossible  to  detect  with  the  FEM  program  in  use.  The  normal  stress, 

Figure  4(b),  is  very  small  for  three-quarters  of  the  width  of  the  line,  reaches  a 
maximum  near  the  edge  and  returns  to  zero,  except  for  at  the  interface  where 
there  is  an  apparent  singularity  as  with  the  other  components.  The  longitudinal 
stress,  Figure  4(d),  is  fairly  constant  throughout  the  width  and  height  of  the  line, 
except  for  the  sharp  rise  near  the  edge  at  the  interface,  and  its  value  is  very 
near  the  biaxial  normalizing  stress. 


Figure  4(a)  Normalized  transverse  stress  at  various  y/h  values  for  an  aspect  ratio  of 
w/h=6.  The  normalization  factor  in  all  cases  is  [E/(1  -v)]AaAT.  All  of  these 


stresses  are  calculated  under  plane  strain  conditions. 


1  2  x/h  3 

Normalized  normal  stress  at  various  y/h  values  for  an  aspect  ratio  of 
w/h=6.  The  normalization  factor  in  all  cases  is  [E/(1  -v)]AccAT.  All  of  these 
stresses  are  calculated  under  plane  strain  conditions. 
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Normalized  Shear  Stress  (w/h  =  6) 


Figure  4(c)  Normalized  shear  stress  at  various  y/h  values  for  an  aspect  ratio  of 

w/h=6.  The  normalization  factor  in  all  cases  is  [E/(1  -v)]AaAT.  All  of  these 
stresses  are  calculated  under  plane  strain  conditions. 


Figure  4(d)  Normalized  longitudinal  stress  at  various  y/h  values  for  an  aspect  ratio  of 
w/h=6.  The  normalization  factor  in  all  cases  is  [E/(1  -v)]AaAT.  All  of  these 
stresses  are  calculated  under  plane  strain  conditions. 
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Figure  5  shows  a  comparison  of  the  transverse,  normal  and  shear  stresses 
at  y/h  =  0.5  for  the  five  different  aspect  ratios.  The  transverse  stress  clearly  gets 
closer  to  the  normalizing  value  as  the  aspect  ratio  is  increased;  this  is  to  be 
expected,  as  a  wider  line  more  closely  approximates  a  film  than  does  a  narrow 
line.  Also,  it  can  be  seen  that  the  percentage  of  the  line  across  which  the  stress 
is  approximately  constant  increases  with  line  width,  being  around  25%  for  the 
w/h  =  2  line  and  closer  to  50%  for  the  w/h  =  10  line.  It  is  evident  that  for  the 
normal  stress  the  influence  of  the  edge  extends  about  one  line  height  into  the 
line,  so  that  for  the  w/h  =  2  line,  the  normal  stress  varies  greatly  across  the 
whole  line,  whereas  for  the  w/h  =  10  aspect  ratio  line  the  stress  is  approximately 
constant  and  close  to  zero  between  x/h  =  0  and  x/h  =  4  ,  and  there  is  only 
significant  variation  between  x/h  =  4  and  x/h  =  5.  Also  note  that  the  maximum 
value  of  the  normal  stress  does  not  change  for  aspect  ratios  greater  than  w/h  = 

4.  The  shear  stress  variations  extend  farther  into  the  lines  than  do  the  normal 
stress  variations,  but  again  the  maximum  value  of  the  stress  does  not  change 
much  for  the  wider  lines.  The  negative  values  of  shear  stress  are  expected, 
since  a  compressive  thermal  strain  is  imposed. 

5.  Conclusions 

The  finite  element  method  has  been  shown  to  be  an  accurate  way  of 
calculating  stresses  in  thin  film  structures.  Furthermore,  while  a  fairly  simple 
geometry  was  used  in  this  paper,  it  would  not  be  difficult  to  extend  the  method  to 
complicated  structures  and  to  include  passivations  and  compliant  substrates. 
Numerical  solutions  to  problems  of  this  sort  would  be  extremely  cumbersome. 

It  has  been  shown  that  as  line  widths  decrease,  thermal  stresses  vary 
across  larger  percentages  of  the  line,  so  that  for  a  w/h  =  2  line  the  stress  cannot 
be  taken  as  constant  across  any  of  the  line.  The  current  frontiers  of  integrated 
circuit  technology  require  smaller  and  smaller  line  widths,  so  it  is  ever  more 
important  to  learn  about  the  stresses  in  these  lines.  The  current  calculations 
should  provide  useful  information  about  the  distribution  of  these  stresses. 


0  1  2  x/h  3  4  5 


Figure  5(b).  Comparisons  of  normal  stresses  at  y/h  =  0.5  for  five  different  aspect 
ratios. 
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Figure  5(a).  Comparisons  of  shear  stresses  at  y/h  =  0.5  for  five  different  aspect 
ratios. 
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